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Abstract 

We give a canonical form for a complex matrix, whose square is nor- 
mal, under transformations of unitary similarity as well as a canonical 
form for a real matrix, whose square is normal, under transformations 
of orthogonal similarity. 

Keywords: Canonical matrices; Normal matrices; Unitary similar- 
ity; Orthogonal similarity. 

For a complex matrix A such that A 2 is normal, we have previously 
used a regularization algorithm and the theory of cosquares to construct two 
canonical forms under transformations of unitary similarity A i— > S~ 1 AS (S 
is a unitary complex matrix) [2\. We now obtain these canonical forms by 
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employing ideas of Littlewood's algorithm [31 [5] . For a real matrix B whose 
square is normal, we also give a canonical form under transformations of 
orthogonal similarity B \— > R~ X BR (R is an orthogonal real matrix). 

Our results provide canonical matrices of linear operators A : U — > U such 
that A 2 is a normal operator and U is a unitary or Euclidean space since 
changes of the basis transform the matrix of A by unitary or, respectively, 
orthogonal similarity. 

The problem of classifying linear operators A : U — > U such that A 3 = 
and U is a unitary space contains the problem of classifying arbitrary linear 
operators on unitary spaces (see [5J p. 45] or |2J). Thus, the condition "A 3 
is normal" and even the condition U A 3 = 0" do not simplify the problem of 
classifying linear operators A on a unitary space. 



1 Squared normal complex matrices under 
unitary similarity 

Theorem 1 ([2]). Let A be a square complex matrix such that A 2 is normal. 
Then A is unitarily similar to 



(a) a direct sum of blocks, each of which is 



or 



[i r 
-fi 



A, /x e C, ^ arg(/x) < 7r, 
r 6 R, r > 0; 



and 



(b) to a direct sum of blocks, each of which is 



or t 



1 
v 



\,uEC, \v\ < 1, 
r el, r > 0. 



(2) 



These direct sums are uniquely determined by A, up to permutation of sum- 
mands. Conversely, if A is unitarily similar to a direct sum of blocks of the 
form ([1]) or (j2J), then A 2 is normal. 

Proof, (a) Let A be a squared normal complex matrix. Let Ai, . . . , At be all 
the distinct eigenvalues of A ordered such that 

< arg(Ai) < 7T and - A* e {A x , . . . , Aj =>- X i+1 = —A, (3) 
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for every nonzero eigenvalue Aj. Schur's unitary triangularization theorem 
[H Theorem 2.3.1] ensures that A is unitarily similar to a matrix of the form 



T 



'A Al Ti 2 



it 



A 



A, 







T-t-\,t 
Aa ( . 



in which every A\. is an x rij upper triangular matrix of the form 

"A,; 







A; 



If is an eigenvalue of A and the corresponding block A = [a^] is nonzero, 
then the sets of indices {i \ ajj ^ 0} and {j \ ^ 0} are disjoint because A 
is upper triangular and A 2 , = 0. We reduce A by permutation similarity 
transformations to the form 



Ar 



* 





with square diagonal blocks. 



Since A 2 is normal, T 2 is normal, too. But T 2 is upper triangular, hence 
T = X^I ni © \2In2 ffi ■ ■ ■ ffi A + In 



H ± n t - 



This implies that K.\ i = \I ni if Aj 7^ 0. 

If Tij 7^ 0, then A 2 = A 2 since T commutes with T 2 . By (j3J) we have 
j = i + 1, < arg(A,) < tt, and A i+ i = -A*. 

Thus, T is a direct sum of matrices of two types: 



XI and T„ 



-fil 



(0 < arg(/i) < tt, F„ ^ 0). 



If = WE^V is a singular value decomposition, then T M is unitarily similar 

to 



s-%s 



-/if 



S :-- 



u 
v- 1 
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Therefore, A is unitarily similar to a direct sum of matrices of the form (pQ). 
Let us prove that this sum is uniquely determined by A, up to permutation of 
summands. We give a direct proof, though we could use Theorem 3.1 in [5], 
which states that every system of linear mappings of unitary spaces uniquely 
decomposes into a direct sum of indecomposable systems, up to isomorphism 
of summands. 

Let T and T' be two matrices that are direct sums of blocks of the form 
(pQ). Let A be an eigenvalue of T. Grouping together the summands with the 
eigenvalues A and —A, we obtain 

T = T 1 ®T 2 , T' = T[ © Tg, 

in which Ti and T[ are direct sums of blocks whose eigenvalues are A or 
—A; T 2 and T 2 ' are direct sums of blocks that have no eigenvalues A and —A. 
Suppose T and T' are unitarily similar. We claim that their decompositions 
into direct sums coincide up to permutation of summands. It suffices to 
verify that the decompositions of Ti and T[ coincide up to permutation of 
summands. Let S be a unitary matrix such that TS = ST'. Since T\ and 
T[ have no common eigenvalues with T 2 and T 2 ', the matrix S has the form 
S = Si © 5 2 and 

715-1 = S X T[. (4) 

Thus, T\ is unitarily similar to T[. 

If Ti is a direct sum of 1 x 1 blocks [A] and [—A], then the decompositions 
of Ti and T[ coincide up to permutation of summands. Suppose Ti has a 
2-by-2 direct summand of the form ([1]). Reduce Ti and T{ by permutation 
similarity transformations to the form 



Ti 



XL 



-XL 



XIpi 



-xi q/ 



in which 



E = diag(ri, . . . ,r k ) © 0, 
are real matrices such that 

n ^ . . . ^ r k > 0, 



E' = diag(r^, 



r[ ^ • • • ^ r v > 0, 



and if A = then E and E' have no zero columns (hence their columns 
are linearly independent; we gather the zero columns and rows in the first 
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vertical strip and the first horizontal strip of T\ and T[). Then (J3J) implies 
that S\ = U © V, in which U is p X p. Thus, SV^ = UT,' and the uniqueness 
of a singular value decomposition ensures that S = £'. 

We have proved that all direct sums of matrices of the form ([1]) are canon- 
ical under unitary similarity. 

(b) Let us prove that all direct sums of matrices of the form ([2]) are 
canonical, too. It suffices to verify that the mapping 



T 



"0 


1 




\fy 1 — \v\ 









. ~V" _ 


V 



(in which ^ arg(-y/z/) < it; that is, \Jv is the principal square root of v) is 
a bijection of the set of matrices of the form 



M U , T := r 
onto the set of matrices of the form 



1 
v 



N, 



v e C, \v\ < l, 

tGI, r > 0, 



\i 6 C, ^ arg(yu) < 7r, 
rel, r > 0, 



(5) 



and each M V<T is unitarily similar to /(M„ iT ). 

First we prove that / is a bijection. Fix N^ r and verify that it has exactly 
one preimage M v<r . The equality j{M v<T ) = r is valid if and only if 



if and only if 



T\JV 

r 



fx, 



r(l — \u\) = r 



v 



(6) 
(7) 



Write the complex numbers fi/r and v in polar form: fi/r = pe tip and v 
Xe 1 ^. The first equality in ((7|) is valid if and only if 



P 



y/X 



t 
2' 



(8) 



There is a unique real x i n the interval (0, 1) that satisfies the first equality 
in (ED because the real function 



x 



x 
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steadily increases from to +00 on the interval [0, 1). Thus, there is a unique 
v that satisfies the first equality in ([7]). We find r from the second equality 
and obtain the required preimage M V ^ T of N^ r . 

Each matrix pair (M U>T , N^ r ) with N^ r = /(M„ iT ) is completely deter- 
mined by the parameters fi G C and r G R. These parameters satisfy 
^ arg(/i) < 7r and r > since by ([H]) they determine the remaining param- 
eters v and r: 



v 



r(l - \u\) = r 



T 



Using these equalities, we obtain 

/i(r -r) 



n 2 n 



V[LT 



and so 



S:~- 



■sJt 2 + /i/i L ^ 



T 



Therefore, M^ T is unitarily similar to N, 



□ 



Littlewood's algorithm [31 El [5] transforms each square complex matrix 
A into a matrix A can that is unitarily similar to A. Two square matrices 
A and S are unitarily similar if and only if they are transformed into the 
same matrix A can = B can . Thus, A can is a canonical form of A under unitary 
similarity. The structure of A cail is studied in [HE]. If A is squared normal, 
then A can is permutationally similar to a direct sum of blocks of the form 



2 Squared normal real matrices under or- 
thogonal similarity 

The realification of an m x n complex matrix M is the 2m x 2n real matrix 
M R obtained by replacing every entry a + hi of M by the 2x2 block 

a —b 

b a 

The real Jordan form of A e M nxn can be obtained from the canonical 
Jordan form of A by replacing all pairs of complex conjugate Jordan blocks 

J n (a + bi) © J n (a — bi), b > 0, 
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by J n {a + bi) R , see [TJ Theorem 3.4.5]. A real canonical form of a normal 
matrix A G IR nxn under similarity can be obtained from its diagonal canonical 
form by replacing all pairs of complex conjugate diagonal entries 



[a + bi] © [a -hi], 



b>0, 



by [a + 6i] R , see [TJ Theorem 2.5.8]. In the following theorem we show that a 
real canonical form of a squared normal matrix A e R riXn under orthogonal 
similarity can be obtained in the same way from the canonical form of A 
given in Theorem QJb) (and from any other canonical form of A under unitary 
similarity) . 

Theorem 2. Let A be a square real matrix such that A 2 is normal. Then A 
is orthogonally similar to a direct sum of real blocks of the form 



T 



H<1, r>0), 



(9) 



and 



a —b 
b a 



(6>0), 
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1 

c -d 

d c 



(d > 0, c 2 + rf 2 < 1, r > 0). (10) 



77ms direct sum is uniquely determined by A, up to permutation of summands. 
It can be obtained from the canonical form of A under unitary similarity that 
is a direct sum of matrices of the form ([2]) by replacing all pairs of summands 



[a + bi] © [a - bi] (b > 0), r 








1 




1" 


r 


c + di 





©T 


c — di 



(d>0) 



by the corresponding matrices (11 01) . Conversely, if A is unitarily similar to 
a direct sum of blocks of the form Q and (jlOO . then A 2 is normal. 

Proof. Let us prove that a complex matrix M is squared normal if and only 
if its realification M R is squared normal. If M is represented in the form 
M = A+Bi with A and B over R, then its realification M R is permutationally 
similar to 

M * == \ A B ~A 
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Since 



A + Bi 




~I 


il ' 




I 


il ' 




'A -B 


A — Bi 




/ 


-il 




I 


-il 




B A 



we have 

Mr = S~\M © M)S = S*(M © M)S, 

in which 



is unitary. 

Thus, M R is squared normal if and only if Mr is squared normal (i.e., 
(M 2 )*M 2 = M 2 (M 2 )*) if and only if 

(M 2 )*M 2 = M 2 (M 2 )* and (M 2 )*M 2 = M 2 (M 2 )* 

if and only if M is squared normal. 

A canonical form of real matrices under orthogonal similarity can be ob- 
tained from any canonical form of complex matrices under unitary similarity 
as follows. We say that a square complex matrix A is decomposable if it is 
unitarily similar to a direct sum of square matrices of smaller size. Let S be 
any set of indecomposable canonical complex matrices under unitary similar- 
ity (for example, the set of indecomposable matrices on which Littlewood's 
algorithm acts identically). Each matrix M e S that is unitarily similar to 
a real matrix R, we replace by R. Each pair {M, N} C S in which M is not 
unitarily similar to a real matrix and N is unitarily similar to the complex 
conjugate matrix M, we replace by M R or by N R (N = M if M is unitarily 
similar to M). Denote the set obtained by <Sr. Theorem 4.1 in [5] about sys- 
tems of linear mappings on unitary and Euclidean spaces ensures that each 
real matrix A is orthogonally similar to a direct sum of matrices from 5r and 
that this sum is determined by A uniquely up to permutation of summands. 

Let us remove from the set S all matrices that are not squared normal 
and construct the set <Sr as above. Then <Sr consists of squared normal 
real matrices and each squared normal real matrix is orthogonally similar 
to a direct sum of matrices from «Sr, which is determined uniquely up to 
permutation of summands. 

Let us prove that if S is the set of matrices (T5]), then 5r is the set of 
matrices (Q and ([TO]) . The real matrices of the form (T5]) give (Q. Each pair 
{[a + bi], [a — bi]} C S with b > gives [a + bi] R , which is the first matrix in 



I il 
I -il 
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0. Each pair {M C+ ^ T , M c _^ r } C S of matrices of the form fl5]) with d > 
and c 2 + d 2 < 1 gives M^. dj T , which is the second matrix in (1101) . □ 
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